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ABSTRACT 

We study the decay process of an unstable quantum system, especially the de- 
viation from the exponential decay law. We show that the exponential period no 
longer exists in the case of the s-wave decay with small Q value, where the Q value 
is the difference between the energy of the initially prepared state and the minimum 
energy of the continuous eigenstates in the system. We also derive the quantitative 
condition that this kind of decay process takes place and discuss what kind of system 
is suitable to observe the decay. 



1 



I INTRODUCTION 



Since the period of the classic works by Dirac [I] and Weisskopf and Wigner [2] , it 
has been a problem how to describe the decay process of an unstable state following 
the principles of quantum mechanics. As is well known, the survival probability of 
the initial state P(t), which concerns the decay of the quantum state, is frequently 
described by the exponential decay law P(t) = e~ rt . However, it is also known that 
the decay process does not obey the exponential law precisely, so it has always been 
a question of how the deviation from the exponential decay law occurs, especially at 
the late and early times of a decay process 

Theorists are always motivated to work on this old problem when high-resolution 
experiments, which are accomplished by a new technology, are performed to detect 
the deviation of the decay law from the exponential [3] - In addition, recent 
several experiments have reported the measurement-induced suppression in quantum 
systems at the early stage of decay, which may be a result of the Quantum Zeno Effect 
(QZE) 0. 

As mentioned above, the deviation from the exponential law at late and early 
times is often discussed [HIE]- At very late times, the survival probability P(t) must 
decrease more slowly than the exponential and exhibits the inverse power law of 
time P(t) ~ t~ a , where a is positive and depends on the property of the unstable 
system. At early times, the survival probability decreases following a Gaussian law 
(the square of time t 2 ), which appears inevitably in a quantum process (and causes 
the QZE). Thus the decay of the unstable state proceeds through three stages in 
general. The initial stage is characterized by a Gaussian law, the intermediate stage 
by an exponential law, and the final stage by an inverse power law. 

In this paper we focus on a different mechanism of deviation from the exponential 
law. Such a decay process occurs in the case of small-Q-value s-wave decay (SQS 
decay). Here the Q value is defined by the difference between the energy of the 
initially prepared state (denoted by E ) and the minimum energy of the continuous 
eigenstates (denoted by E th ) in the system. The small- Q- value decay has been 
discussed in some papers [§J EH]- The point we would like to emphasize here is as 
follows: in the case of the SQS decay, we can observe not only the enhancement of 
the QZE, but also no exponential period. This means that the deviation from the 
exponential law can be observed easily if the SQS decay system is prepared. We also 
derive the quantitative condition that such a decay takes place. 
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Figure 1: The shape of the potential. 



This paper is organized as follows. In Sec. II, we show an example of the SQS 
decay by using the one-dimensional tunneling system with a box-type potential. 
The general description for an unstable system is formulated in Sec. III. Using this 
formalism, we derive the quantitative condition that the SQS decay occurs in Sec. IV. 
In Sec.V, we summarize our results and discuss what kind of system exhibits such a 
decay process. 



II AN EXAMPLE OF SQS DECAY IN TUNNELING PHENOMENA 

Before going into the general discussion, we show the example of SQS decay in 
the tunneling phenomena. In this section we discuss the one- dimensional tunneling 
problem because only the radial part of the wave function is relevant to the s-wave 
tunneling even in a three-dimensional system. 

Let us consider the decay process through the one- dimensional box-type potential 
depicted in Fig.l. Parameters for characterizing the system are also shown in the 
figure. Here we assume U is not so large that there is no bound state. The goal 
here is to calculate the survival probability of the prepared state and show that the 
SQS decay is realized in this system. 

The survival probability P(t) is defined by the nondecay amplitude a(t) as P(t) = 
|a(t)| 2 . a(t) is given by 

a(t) ee (0|e^|0) , (1) 
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where |0) is the initially prepared state and H is the Hamiltonian of the system. 
Please note that we use the units H = 1 in this paper. The nondecay amplitude can 
be expanded over the energy eigenstates \E) as 

a(t)= / dE p(E) e~ iEt , p(E) = \(E\0}\ 2 . (2) 
Jo 

The function p(E) is called the spectral function, in which all information about the 
decay process is included. The energy eigenstate \E) can be obtained analytically 
in this system as detailed in Appendix A. Here we take the initial state |0) as the 
ground state in the well for the infinitely height barrier, 

il>i(x) = (x\0) = j| sin {^)0{x + a)0(-x) . (3) 

The energy expectation value of this state is E Q = ir 2 /(2ma 2 ) — Uq and the Q value 
is given by Q = Eq — E t h = tt 2 / (2ma 2 ) — Uq because the spectrum of the continuum 
energy eigenstate starts from the zero energy. Using this initial wave function, the 
spectral function of the system is obtained analytically after some calculations, and 
given by 

1 27rgsin 2 r 
P[ ' ~ 2ma 2 a{E) (r 2 - vr 2 ) 2 ' U 

where 

a(E) = q 2 + Gq cos 2 r sin 2 r H ( - cosr • sinh(sw) + sinr • cosh(s-u) I , (5) 



u u 

for the case that the energy is smaller than the potential barrier E < U/b, while 

G G (t \ 2 

a(E) = q 2 + G cos 2 r sin 2 rH (- cosr • sin(sw) + sinr • cos(sw) ) . (6) 

u u \s J 

for the other case E > U/b. Here we use the dimensionless quantities to write down 



the spectral function q = ^2mEa, r = J2m(E + U )a, s = J2m(U/b — E)a, and 



s = ^2m(E — U/b)a. Furthermore, the variables that characterize the potential are 
also given by the dimensionless ones G = 2maU, Gq = 2ma 2 [/ , and u — b/a. 
For investigating the SQS decay, we define the decay rate T(t) by 

r(i) = -| In P(t). (7) 

This quantity is more convenient rather than the survival probability itself because 
this rate is constant while the decay process is governed by the exponential law. We 
performed the integral in Eq.(J2J) numerically using the spectral function Eq.(@J), and 
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Figure 2: Examples of decay. Here we set G = 20, G = (or Q = 8.97365(2ma 2 )" 1 ), 
u = 10~ 4 on the left and G = 20, G = 8.957335 (or Q = 6.55445 x lO^ma 2 )- 1 ), 
u = 10 -4 on the right. Each quantity in the figure is averaged over a short time interval. 
The example on the right corresponds to the SQS decay. 

calculated the decay rate T(t). The results are shown in Fig.2. We studied two cases 
of Go, that is, Go = and 8.957335 which correspond to the cases Q = 8.97365 and 
6.55445 x 10~ 4 [the units of Q are (2mo 2 ) -1 ], with fixed G and U. 1 As you see, the 
exponential decay law is observed in the case that the Q value is not very small. 
On the other hand, if the Q value is small enough the exponential period no longer 
exists even at the time t ~ 500(2ma 2 ) when P(t) decreases to the order of e -1 . This 
is nothing but an example of the SQS decay. 

In the following sections we investigate the SQS decay by a general description 
and what kind of situation is necessary. We also derive the quantitative condition 
for the SQS decay to take place. 

Ill GENERAL DESCRIPTION FOR UNSTABLE STATE 

In this section we explain the general formalism of unstable state decay. Using 
this formalism the quantitative condition for the SQS decay is derived in the next 
section. 

For calculating the nondecay amplitude we use an exact integro-differential equa- 

1 Rigorously speaking, there is a quantum correction to the energy level and Q should be given 
by eq.(26), where Q is indeed defined in this sense. However, this correction is assumed to be small 
and Q is set to be Eq — Eth in other lines or equations. 
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tion by using the technique in Refs. [TTJ We introduce the projector onto the 
initial unstable state, V = |0)(0|, and decompose the Hamiltonian H as 

H = H Q + V, 

H = VHV + (1 - V)H(1 - V) , V = H-H . (8) 

Define the energy eigenstates \a) of Hq; Ho\a) = E a \a). It is then easy to confirm that 
the decay interaction V operates only between the initial state |0) and its orthogonal 
complement |n); Vq u = V* ^ 0, V mn = V 00 = 0. Here we use the intermediate roman 
letters such as m, n to denote eigenstates projected by 1 — V. The interaction V 
depends on the initially prepared state. Although this formalism looks odd, we can 
execute an exact analysis like a time development of the unstable state using this 
tool. 

We work in the interaction picture and expand the state at a finite time t, using 
the basis of the eigenstate of H ; \ip)i = e lHot \ip) = S a c a (t)|a) . We thus write the 
time evolution equation for the coefficient c a (t): 

ico = E V °n e-^~ E ^ c n , ic m = V m0 e ^- E °)* c . (9) 

n 

Here Eq is the energy of the initial unstable state; E$ = (0\H\0). The nondecay 
amplitude is related to this coefficient by (0|e~* H *|0) = e~ lEot co(t). From the above 
equations, a closed form of the equation for the nondecay amplitude c (t) = a{t) 
then follows: 

a{t) = - f dt P{t - t') a{t') , (10) 



f3{t-t') = (0|V7(t)V}(t')|0) = / dE a{E) e -^-E )(t-t>) ? (n) 
a{E) = J2 5 (E-E m )\V 0m \ 2 . (12) 



in 



Here Vj(t) = e lHot V e~ lHot is the decay interaction written in the interaction picture 
and the function cr(E) characterizes the interaction between the unstable state |0) 
and the other states \m). The initial condition c m (0) = is used to derive the 
equation for a(t), and E t h is the threshold for the state \m). 

The standard technique to solve this type of integro-differential equation (jlOj) 
is the one that utilizes the Laplace transform, and we finally obtain the nondecay 
amplitude in the form 

a(t) = r ^F(E + zO + ) e-^-Eo)' , (13) 

J-oo ZTTl 

1 r°° a(E) 

F(z) = -— , G(z) = / dE^- . (14) 

y ' -z + E - G{z) ' 1 ; J Eth E-z K ' 
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The initial condition a(0) = 1 is imposed in this derivation. 

The analytic property of the function F(z) is evident; this function is analytic 
except on the branch cut which runs from the threshold value Efh to positive infinity 
on the real axis (Fig.3). As is well known, if the Riemann surface is considered by 
analytic continuation through the branch cut (and regarding the original complex z 
plane as the first Riemann sheet), there is a pole on the second Riemann sheet near 
and below the real axis if the decay interaction is weak enough. The pole location 
z p is determined by 

r°° o~(E) 

z p -E + G n {z p ) = z p -E + / dE v _ y + 2ma(z p ) = , (15) 

where the analytic function G(z), and hence F(z), is extended into the second sheet 
by Gu(E — i0 + ) = Gj(E + i0 + ) through the branch cut. The real function &(E), 
which was originally defined for real E > E t h, is also extended to the function o~(z) 
defined on the complex plane by analytic continuation. In addition to this pole 
there may be some singularities on the second Riemann sheet, but we ignore the 
effects of such singularities in the following discussion. This approximation is valid 
as discussed in Ref. |§] because these singularities do not affect the decay phenomena 
except for the early stage. 

Using the discontinuity of the analytic function F(z) across the branch cut on 
the first Riemann sheet, 

F(E + i0 + ) - F(E - i0 + ) = 2ma(E) \F(E + i0 + )\ 2 , (16) 

which is called the "elastic" unitarity relation, we can deform the contour of integra- 
tion on the real axis in Eq. (jl3|) into the sum of two contours, one around the pole 
as shown by C p and the other along C t h in Fig.3, 

a(t)=(f +[ ) ^F(z) e-**-** . (17) 

We consider the case that the unstable state initially prepared is a metastable 
state, which means the decay interaction is weak in comparison with the typical 
energy of the system, E sys , induced from the oscillation in the well: cr E sys ~ E . 
Then the pole location on the second Riemann sheet can be obtained approximately 

as 
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Figure 3: The analytic structure of the complex z plane and the contours of the integrals. 
The contours and the pole shown by the broken line are on the second Riemann sheet. 



where p. v. means the principal value of the integration. The C p integration in 
Eq. (jl7|) can be performed without difficulty by the residue theorem, and with the 
aid of the approximation above, this becomes 

f —F(z) e-^-^t ~ e m{E )t e -Ka{E )t _ 
2ni y ' 

This integration gives an 0(1) contribution to the nondecay amplitude a(t). On 
the other hand, the integration along Cth is of O(a/Eo), which gives only a small 
contribution. This is because the integration can be approximated by 

f rj 2~ roo 

/ —F(z)e-^- E ^ = -ie^- E ^ dya(E th -ty)\F(E th -ty)\ 2 e-y\ 
JCth 2m Jo 

~ -ie^-^)*|F(^)| 2 / dy a(E th - iy) e* , (20) 

Jo 

for sufficiently large t, and the factor \F(E th )\ 2 usually takes the value (E — E th )~ 2 . 
The dominance of the C p integration in Eq. (jl7)l leads to the exponential decay 

law 

P(t) = \a(t)\ 2 ~ exp(-r p t) , T p = 2na(E ) , (21) 

and this coincides with the familiar golden rule of perturbation theory. From this 
investigation, we know that the perturbative calculation gives satisfactory results in 
many cases. We would, however, like to elucidate the time evolution in finer detail, 
and study the conditions for breaking the exponential decay law, especially the SQS 
decay, in the next section. 
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IV THE CONDITION FOR THE SQS DECAY LAW 



We now derive the quantitative condition of SQS decay. We also discuss the 
situations in general that exhibit deviations from the exponential decay law. 

The C p integration around the pole z p in Eq. (|17jl always yields exponential time 
dependence, so nonexponential decay is realized when the C t h integration contributes 
to the nondecay amplitude a(t) by the same order as the C p integration. 

From the evaluation of the contour integrations in Eq. (jl7|) . we can classify the 
nonexponential decays that satisfy the condition mentioned above into three cases. 
In the following we describe them in detail. 

The first one concerns the short time behavior and is known as the QZE. At 
early times (t < E~y S ), the approximation used in Eq. (j20j) is no longer valid and 
the high-frequency component of cr(E) becomes important. From the definition of 
the survival probability, we naively expect that the short time behavior exhibits a 
deviation from the exponential law, which is in the form of 

|(0|e- im |0)| 2 ~ 1 - t 2 ((0|# 2 |0> - (0|#|0> 2 ) . (22) 

Thus quantum mechanics appears to predict a quadratic form of deviation in the 
t -> limit. 

The second one relates to the long time behavior. At late times t ^> 1/T P the 
C p integration is exponentially suppressed, while the C t h integration is not strongly 
suppressed because the behavior of <r(E) near the threshold is expressed by cr(E) = 
c(E — E t h) a , which leads to power law behavior of the Cth integration, 

J ^n*) * -ie i(E °- E ^\F(E th ) f- T(a ^l^ 12 ■ (23) 

Here T(z) is Euler's gamma function. Therefore the contribution from the Cth in- 
tegration exceeds the one from the C p integration at late times, and the decay law 
changes from exponential to an inverse power law. 

The small Q (= Eo — E t h) value case is the last one on which we shall mainly focus 
in this paper, that is, the SQS decay. This case can be understood by investigation 
of the prefactor \F(E th )\ 2 in Eq.® in detail: 

\F(Et h )\ 2 = [Q-U(E th )r 2 . (24) 

Since the function U(E t h) is of 0(a), this factor gives Q~ 2 when the Q value is not 
very small. However, if the Q value is of the same order as a, the factor \F(E t h)\ 2 
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Figure 4: The typical shape of the spectral function p(E). The SQS decay is realized 
when the Q value is smaller than the width of the peak. 

becomes large and the contribution from the C t h integration becomes comparable 
with the one from the C p integration. In this case decay that does not include the 
exponential period at all can be realized, and this situation never occurs in other 
cases described above. 

We now derive the quantitative condition that the SQS decay takes place. By 
virtue of the analytic property of F(z), we can express the amplitude a(t) in the 
convenient form 

a(t) = Q d E p(E) e-^ , , (B) = __^__. (28) 

The function p(E) is the spectral function as mentioned in Sec. II. The schematic 
shape of p(E) is shown in Fig.4. This function takes a real positive value when E 
is larger than the threshold E t h- p{E) drops quickly in the limit E — > oo because 
this function must satisfy the normalization condition fjf dEp(E) = 1, which is 
equivalent to the initial condition a(0) = 1. In addition, p(E) has a peak around 
E ~ Eq. (More precisely, the exact location and the width of the peak are determined 
by the real and the imaginary parts of z p in Eq. (|15jl .) Therefore the dominance of the 
C p integration is equivalent to the Breit-Wigner form of p(E) or the limit E t h — * — oo, 
as you can find from the shape of p(E). 

The absence of the exponential decay law due to the small Q value then occurs 
when Q is smaller than the width of the peak because the Breit-Wigner shape of 
p(E) does not hold in this situation. Thus the condition for the SQS decay is given 
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by 

h "'"'>!, Q = ftz p -E th . (26) 



Q 

This is the final result of this section, the quantitative condition for the SQS decay 
to take place. In Ref. [Hj, the authors derived the opposite condition to Eq.(J2H} 
from the viewpoint that nearly exponential decay takes place. The condition here is 
a necessary and sufficient condition, so our condition is consistent with theirs. 

To investigate the condition Eq. (|26|) for a given system, this inequality is not very 
useful because the calculation for the exact location of z p from Eq.(|15|) is a tedious 
one. This situation is somewhat improved if we use the approximation Eq. (jl8|) ; then 
the condition turns to 

Tp > 2 . (27) 



Eq — E th 

This is more convenient for practical use. 

As an example, we show the contour plot of ^sz p /Q for the model discussed in 
Sec. II, which is depicted in Fig.5. The contour plot is drawn on the (G, Q) plane 
with fixed u. The crosses in the figure correspond to the values used in Fig. 2, and 
in the shaded region the parameters satisfy the condition that the SQS decay takes 
place. 

V DISCUSSION AND SUMMARY 

We now come to the stage of discussing what kind of system is necessary for 
the SQS decay. The important point we must pay attention to for discussing the 
condition Eq. (}2T)j) is the Q dependence of Qz p , the imaginary part of the pole location. 
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The Q dependence is determined by the function o~(E ) as shown in Eq. lfTHj) . When 
the Q value is small enough, the dependence is determined by the threshold behavior 
of a(E ), which originates from the spectral function p(E ) as shown in Eq.()25J). 
Therefore the threshold behavior of the spectral function, taken to be p{E ) = c(E Q — 
Eth) a = cQ a , is the key quantity in this problem. Here the coefficient c is a constant 
dependent on the system. The SQS decay is expected when the power of p{E ) near 
the threshold, a, is smaller than 1. 

As is well known, the threshold behavior of the spectral function is determined by 
the quantum number of the orbital angular momentum in a scattering process (e.g., 
a particle decay or a radioactive process), which we denote / jT^j. The threshold 
behavior is given by p(E ) = c(E — E t h) l+1 ^ 2 - Therefore s-wave (I = 0) decay is 
necessary for the SQS decay in these processes, and the SQS decay never occurs via 
higher-/ (/ > 1) processes. 

Let us move on to the decay process through tunneling. For the one-dimensional 
model discussed in Sec. II, the threshold behavior is given by 

= sSo) ' Q = E a - E lh , (28) 

which is the same as in the case of the s-wave decay. This is the very reason that 
the system exhibits the SQS decay when the Q value is small enough. You might 
think that such a threshold behavior is due to a peculiarity of the potential. This 
is, however, not correct. To check this, let us consider a system with the modified 
potential shown in Fig. 6. The spectral function of this system can also be obtained 
analytically using Bessel functions. The threshold behavior again coincides with 
the case of the s-wave decay. The form of the spectral function and the threshold 
behavior are given in Appendix B. 

This result is naturally understood if we consider the model with spherical sym- 
metry in three dimensions. Since the angular momentum is a conserved quantity 
in this case, the s-wave decay process can always be reduced to a problem in a 
one-dimensional system. Therefore the threshold behaviors of all models in one 
dimension with nonsingular potential are the same as the one of s-wave decay. 

We summarize the results of this paper. The decay of the unstable state with 
small Q value in the s-wave process (SQS decay) exhibits the interesting feature 
that there is no exponential period. If we can make practical use of this mechanism, 
we will easily observe the deviation from the exponential law. As mentioned above, 
the SQS decay takes place in a system described by an s-wave process or a one- 
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Figure 6: The modified potential. 

dimensional system. However it is difficult to prepare a setup of small Q value in 
experiments on particle decay or radioactive processes, because the Q values in such 
cases are fixed by nature and we cannot control them. On the other hand, the 
tunneling phenomenon may be hopeful to observe the SQS process because we may 
achieve sufficiently small Q more easily. Unfortunately, experiments to look for the 
SQS decay have not been carried out until now. It is important to discuss an actual 
physical system which realizes the SQS decay, but this issue is beyond the scope of 
this paper and remains as a future problem. We believe that this kind of experiment 
is interesting to observe the nonexponential decay of an unstable quantum system. 
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APPENDIX A 

In this appendix, we calculate the spectral function p(E) = \(E\0)\ 2 of the one- 
dimensional model with a box-type potential that is used in Sec. II. The explicit form 
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of the energy eigenstate \E) is necessary in this computation and it can be obtained 
analytically 

{A sin r(x/a + 1) for (— a < x < 0) , 

B sinh(sa;/a) + C cosh(sx/a) for (0 < x < b), (29) 

Dsm(qx/a) + Ecos(qxfa) for (b<x), 

for E < U/b and 

A sin r(x/a + 1) for (— a < x < 0) , 

= (x\E) = \ Bsin(sx/a) +Ccos(sx/a) for (0 < x < b), (30) 

D sm(qx/a) + E cos(qx/a) for (b<x), 

for E > [7/6. All parameters such as a, 6, g, r, s, and 5 are defined in Sec. II. 
The coefficients A, B, C, D, and E can be determined by the junction condition at 
x = —a, 0, b and the normalization condition of the eigenstates, (E\E f ) = S(E-E'). 
For example, the coefficient A is given by 

\a\ 2 = q i"*-n 

11 2maHa{E) ' 1 ' 

and others can be obtained similarly. The function a(E) is defined in Eqs.(jSJ) and 
(JHJ). The spectral function is defined by the overlap of the energy eigenstate \E) and 
the initially prepared state |0) given in Eq.(JSJ). Using the eigenstate obtained above, 
the spectral function is computed as 



p(E) = \(E\0)f 



dx<j} E (x)*ipi(x) 



2 1 27rg sin 2 r 

2ma 2 a(E) (r 2 — tt 2 ) 2 



(32) 



APPENDIX B 



In this appendix, we write down the explicit form of the spectral function of the 
model whose potential is depicted in Fig. 6 and which is used in the discussion of 
Sec.V. The spectral function of this system can also be obtained analytically using 
Bessel functions, and the computation of the spectral function can be performed in 
the same way as in Appendix A. After some calculations we obtain 

(F] = 1 7rsin 2 (fcia) 

P{ ' 2ma^k z {\e 1 \ 2 + \e 2 \ 2 ){k 2 a 2 -T: 2 ) 2 ' 1 J 
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where k\ = \ / 2mE, k 2 = y2m(hi — E), and k 3 = y2m(E — h 2 ). All parameters 
characterizing the potential such as a,...,d and hi,h 2 are defined in Fig.6. The 
variables e\ and e 2 are given by 



F-\k 3 ) ■ T R (z d ) ■ T^(z c ) ■ Ph ■ F-\-ik 2 ) ■ T L (z b ) ■ T L \z a ) * .(34) 




L and R appearing as the subscripts of the matrix T represent the slopes of the barrier 
at the left and the right sides, L = 2ma 3 hi/(b — a), R = 2ma 3 (hi — h 2 )/(d — c). The 
matrices F(k) and T t (z) are defined by 

m = I .! \), T M . ( „f^ W , ±!> {Z } , I . (35) 



ika -ika ] \tV z B' 1/z {z) t 1/3 ^ 1/3 



z 



The arguments of the matrix, z a , . . . , z d are the rescaled positions of a, ... ,d, which 
are defined by z a = —L l ^ 3 (b — a)E/(ahi), z b = L 1//3 (6 — a) (hi — E)/(ahi), z c = 
R 1 / 3 (d - c)(hi - E)/(a(hi - h 2 )), z d = R 1 ' 3 (d - c)(h 2 - E)/(a(h 1 - h 2 )). The 
function B±ya(z) is defined by the Bessel function J±U3 as 



5 ± i/3(^^%f!(-) 3/2 ) > (36) 



3 

and B'(z) means the derivative with respect to z. 5i and S 2 at the right side in 
Eq. (j34|) are the values of the wave function and its derivative at the origin, which 
are given by Si = sin^a), S 2 = kiacos(kia). Finally, the matrix Ph is given by 

/ P k 2 (c-b) p ~k 2 {c-b) \ 

Ph=\ , , . (37) 

y k 2 ae k2 ^ k 2 ae~ k2 ^ J 

In the rest of the appendix we show that the threshold behavior of this system is 
same as that of the s-wave decay. We expand the spectral function with respect to 
(E — h 2 ) in the vicinity of the threshold. The components of the matrices in Eq. (j34|) . 
T R (z d ), T R 1 (z c ), Ph, ik 2 ), T L (z b ), and T£ (z a ), and the variables Si,S 2 become 

constant at the leading order, and only the matrix F~ 1 (k 3 ) has E dependence such 
that 




-ik^ 
-ik 3 1 



(38) 



Therefore ei and e 2 are proportional to k^ 1 . As a result, the threshold behavior of 
the spectral function p(E) of Eq. (j33j) is of the form p(E) oc k 3 = (E — h 2 ) l l 2 , which 
is same as the result of s-wave decay. 
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